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ABSTRACT 

Aims. The Maunder minimum (MM) of greatly reduced solar activity took place in 1645-1715, but the exact level of sunspot activity 
is uncertain as based, to a large extent, on historical generic statements of the absence of spots on the Sun. Here we aim, using a 
conservative approach, to assess the level and length of solar cycle during the Maunder minimum, on the basis of direct historical 
records by astronomers of that time. 

Methods. A database of the active and inactive days (days with and without recorded sunspots on the solar disc respectively) is 
constructed for three models of different levels of conservatism (loose ML, optimum MO and strict MS models) regarding generic 
no-spot records. We have used the active day fraction to estimate the group sunspot number during the MM. 

Results. A clear cyclic variability is found throughout the MM with peaks at around 1655-1657, 1675, 1684 and 1705, and possibly 
1666, with the active day fraction not exceeding 0.2, 0.3 or 0.4 during the core MM, for the three models. Estimated sunspot numbers 
are found very low in accordance with a grand minimum of solar activity. 

Conclusions. We have found, for the core MM (1650-1700), that: (1) A large fraction of no-spot records, corresponding to the solar 
meridian observations, may be unreliable in the conventional database. (2) The active day fraction remained low (below 0.3-0.4) 
throughout the MM, indicating the low level of sunspot activity. (3) The solar cycle appears clearly during the core MM. (4) The 
length of the solar cycle during the core MM appears 9 ± 1 years, but there is an uncertainty in that. (5) The magnitude of the sunspot 
cycle during MM is assessed to be below 5-10 in sunspot numbers; A hypothesis of the high solar cycles during the MM is not 
confirmed. 

Key words. Sun:activity - Sun:dynamo - Sun:Maunder Minimum 


1. Introduction 

There was a period, in the second part of the 17th century, 
of greatly reduced solarjictivity that was named the Maunder 
Minimum (MM) by l Eddv (19761) . The MM was characterized 
by almost complete absence of reported sunspots on the so¬ 
lar surface although some indications of cyclic activity can 
be notic ed particularly in the geomagnetic and heliospheric 
indices dBeer et al J1 1 998b lllsoskin et all 1200 it ISoon & Yaskelll 
120031) . The reconstruction of solar activity based on the his¬ 
torical records_of_telescopic observations of sunspo ts since 
1610 (lHovt & Schattenll998al:lHovt & Schattenl 1998a -called 
HS98 henceforth) marked a milestone in the study of solar ac¬ 
tivity in the recent past and, especially, for the MM period. The 
Group Sunspot Number (GSN) built by HS98 became the only 
high-resolution (daily) index to study solar activity during the 
MM. 


covered by sunspot data and more than 90% of days have formal 
observation records in the HS98 database. However, it contains 
a large number of generic statements of the absence of sunspots 
during a long period of time. Such records are not strict observa¬ 
tional data but they were interpreted by HS98 as no-spot data. 
Many o f these records corresponded to s olar meridian obser¬ 
vations ( Vaau erol 120071: IClette et al1l2014l) and should be used 
with caution for the reconstruction of solar activity, as shown 
by [Vaquero et al.l (12014! ) who analyzed sunspot records taken 
during systematic solar meridian observations performed at the 
Royal Observatory of the Spanish Navy from 1833 to 1840. 
Moreover, as lCarrasco et al.l (120151) suggested basing on an anal¬ 
ysis of sunspot records by Hevelius in the 17th century, the GSN 
index may be underestimated during the MM due to a large num¬ 
ber of ’’zero” sunspot records taken from solar meridian obser¬ 
vations. In general, astrometric observations of the Sun are not 
always reliable for sunspot counting because of the different aim 
of such observations. For example, there is no information on 
sunspots in the extensive table of astrometric records of the Sun 
made with the meridian line in the San Petronio Basilica from 
1645 to 1735 as publish ed bv iManfredil (1 1 7361) . Nevertheless, 
iHovt & S chatted (1199 8a) adopted solar observations recorded in 
this source as no-spot reports, which is not correct. It has been 
discussed that, while the definition of sunspot numbers and even 
sunspo t groups is not very reliable in the earlier part o f the GSN 
series dciette et al.ll2014: Zolotova & Ponvavinll2015l) . solar ac¬ 
tivity during the MM can be reliably represented by the fraction 


The aim of this work is to elucidate whether the absence 
of the sunspot cyclic activity during the MM was real or an 
artefact caused by a problem of compilation of the database 
of sunspot records. Several studies pointed to possible in¬ 
consistences i n the database used by HS98 especially around 
the MM (e.g., Vaquero & Vazquez! 120091: [V aauero et al.l 1201 ll 


I Va quero & Trigo 2014 ). As an extreme, Bolotova & Ponvavin 

( 20151 ) claimed there was no grand Maunder minimum and that 
sunspot cycles during MM were as high as » 100, which is 
higher than the current cycle # 24. We note that the MM is well 
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of active days dKovaltsov et alJ2004t IVaauero et al.ll2012l l2014t 

IIJsoskinll2013h . 


Despite the overall level of activity, the parameters of the 
solar cyclic variability during MM are also important to know. 
Altho ugh the solar cycle was perceptible in the butterfl y dia¬ 
gram dRibes & Nesme-Ribesll 1 993k IVaauero et al.ll2015h based 
on the observations of sunspot latitudes during the last decades 
of the 17th century, the 11-year so lar cycle is only marginally 


detec table in the sunspot numbers (Waldmeier 1961; Mendoza 


1997b with a weak 22-year cycle dominated (Usoskinetal 


2001). On the other hand, some works based on data of 
high-resolution cosmogenic 14 C measured in tree trunks sug- 


Grand solar minima ( 

Stuiver et al. 1998; Mivahara et al. 

2004; 

Mivahara et al. 2006^ 

Mivahara et al. 2010; Nagava et al. 

2012; 


of the solar cycle was increased to about 14 years during the 
MM, to about 13 years in the beginning of the Sporer Minimum, 
up to 16 years during the 4th Century BC Minimum, and to 12- 
13 years during the late 7th century minima. 

In this work we aim to study variability of solar activity dur¬ 
ing the MM using the statistics of the active days basing on only 
th e most reliable solar obse rvations from the database compiled 
bv lHovt & Sc hat ten! ( 1 998 a ) and to establish an uppermost upper 
(maximum maximorum) limit on that. 


2. Sunspot activity database 

Since quantitative interpretation of many records is uncertain 
for that period, we consider only qualitative indicators of the 
sunspot activity for each day for the period 1637-1715 AD. 
Leaving aside the exact number of reported sunspot group in the 
HS98 catalog, we only consider three possible states for each 
day: 

- no-information or missing days; 

- inactive days when we believe there were reliable observa¬ 
tion of the absence of sunspots; 

- active days when at least one sunspot group was explicitly 
reported by at least one observer; 

We build our database of the active and inactive days for three 
models of different levels of conservatism regarding generic 
no-spot records. For the period 1637-1642, we used exactly 
the records listed by IVaquero et al.l (1201 ll) . For the period 
1643-1715, we used the records from the F1S98 database 
(http ://w w w. ngdc. noaa.gov/stp/ sunspot_numbers/group_ 
sunspot_numbers/alldata.dat) for each observer separately. 
In addition, for the year 1672 several active days were added 
according to observations by N. Bion not included into the 
HS98 database (ICasas et al.l [20061) . While the original HS98 
database contains 26508 daily records for the analyzed period 
1637-1715, our models include much less records because of 
rejecting, with different levels of conservatism, generic state¬ 
ment mostly related to no-spot observations. All these models 
provide an overestimated upper bound of sunspot activity due to 
a possible selection bias towards active days. 

2.1. Loose model (ML) 

This model is similar to that by iKovaltsov et al .i (120041) and ig¬ 
nores all the generic statements (longer than a month) in the 
HS98 database, and considers only explicit statements with ex¬ 
act mentioning dates of observations. This affects such generic 


statements as, e.g., by J. Hevelius for no spots during 1645- 
1651, by J. Picard for 1653-1665, by H. Siveriusfor 1675-1689, 
etc. This model is least conservative and is called ’’loose”. It 
includes 13512 observational days which is nearly half of the 
HS98 database. 


2.2. Optimum model (MO) 

The MO model provides a reasonable balance between strictness 
and data acceptance and is considered as the optimum conserva¬ 
tive. For each year, we considered observations of only those ob¬ 
servers who reported at least one sunspot group at any day of the 
year, which would prove that the observer was ’’active”. In this 
way, generic long-extending reports of ”no-spot” were neglected 
but no-spot records of active observers were considered. The 
MO models is biased towards ’’active” years and produces no 
result for the years without sunspot observations. For example, 
if a year is full of definite ”no-spot” records but does not con¬ 
tain a single sunspot observation reported, such year is marked 
as ”no-information” in this model. Alternatively, if an observer 
was ’’active” during a year, his generic ”no-spot” records for this 
year were considered by the model, so that the total number of 
days Nj in the MO model may exceed that for the ML model 
for some years (see Fig.Q}. This model includes 8089 observa¬ 
tional days for the period analyzed, which is roughly */3 of the 
full HS98 database. 


2.3. Strict model (MS) 

In this ’’strict” model we excluded all the generic statements as in 
the ML models, but additionally we treated other no-spot records 
in a very conservative way, so that we consider as inactive only 
days, when at least two observers independently reported that the 
Sun was spotless and there were no other records of sunspots. If 
at least one observer reported sunspots, the day was considered 
as active. All other days were treated as no-information days. 
This is the most conservative approach, especially in the earlier 
part of the Maunder minimum, when the number of documented 
observers was low and they rarely overlap. This model includes 
5159 daily records or Vs of the full HS98 database. 

For each model we define the number of active Na and the 
total number Nt of the accepted observational days in a year. 
Since the annual data are quite noisy (see below) we also con¬ 
sider triennial intervals. In order to keep the strictness, the MO 
model was still operating with annual periods to identify ’’active” 
observers. The results are shown in Fig. |T| for the three models 
as well as for the formal HS98 database. One can see that, while 
the HS98 database covers the entire period pretty well, the three 
models provide a more conservative estimate of reliable observa¬ 
tions, which is greatly reduced in the earlier Maunder minimum 
but is quite solid towards its end. 

An example of the coverage of the data in the three models 
and the formal HS98 database is shown in Fig. [2] for the year 
1676. Although this year was almost fully covered by data in 
the formal HS98 database, except for a short gap in October, the 
three models considered here include much less inactive days 
while keeping the active days. A small discrepancy in the num¬ 
ber of active days is related to excluded interpolations (as in Dec 
22-24) and confusing values (as in Jun 25 when a sunspot record 
by R. Hook was missed in the formal HS98 series) in the HS98 
database. 
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Fig. 1 . The number of active days N A (grey bars) and total ob¬ 
servational days Nj (curves as defined in the legend for the three 
models and the formal HS98 database) per interval used. Panels 
A and B depict annual and triennial data, respectively. 

We emphasize again that the procedure described above 
serves as an uppermost upper bound because of possible over¬ 
suppressing zero-sunspot records. 



Fig. 2. A map of the days allocation for the year 1676 in the 
original HS98 database (denoted as HS) and the three models 
considered here. Each line represents one month (days of the 
month numbered on the top and bottom) for a model. The empty 
white, blue ”0” and red ”1” cells correspond to no-information, 
no-spot and active days. 


3. Results 

3.1. Active day fraction 

From the collected database of sunspot records, we have esti¬ 
mated the fraction of active days F,\ in each model, as follows 
(cf. iKovaltsov et al.ll2004l) . For each interval, either annual or 
triennial, we have a sample of n daily observations with r ac¬ 
tive days reported. Assuming these observation were taken ran¬ 
domly and independently, one can assess the probability of the 
occurrence of exactly s active days within N days during the 
considered interval (a year or 3 years) using the hypergeometric 
probability distribution: 

_ s'. (N - .y)! n'.(N-n)'. 

P(S> ~ (s - r)\ ( N-s-n + r)! ' (n - r)\N\r\ j 

As the optimum value of s* we consider the median value, viz. 
the value of s which yields P(s*) = 2** p(s) = 0.5. The results 
for annual and triennial time intervals are shown, along with er¬ 
ror bars of a 90% (two-sided) confidence interval, in Figures [3] 
and [4] respectively. We note that triennial data were calculated 
from the original daily values using equation (|7]i and not as an 
average of the annual data. 


3.2. Length of solar cycles 


Although the annual data are quite noisy, the triennial ones (see 
TableQ} clearly show a decadal periodicity during the MM. For 
example, Figure[4]suggests maxima of solar cycles around 1639, 
1655-1657, 1675, 1684 and 1705 in all the models. There is 
also an indication of a cycle maximum around 1666 in the MO 
models, but the statistics is low with a single observation for the 
3-year interval. Periods around 1648 and 1693 are poorly known 
with data gaps in the MO model. 

There are four solar activity maxima in the core MM, be¬ 
tween maxima ca. 1657 and 1684. This leads to an estimate 
of the average solar cycle length (max-to-max) during the core 
MM as 9 + 1 years. However, our view of the cyclic evolution 
of sunspot activity during MM is uncertain because of the un¬ 
clear situation around 1648, 1666, and 1693. If we assume two 
hypothetical missing solar maxima during these periods, as e.g., 
Waldmeieri (196lh proposed a cycle maximum in 1649, while 
Usoskin et alj ( 200lb suggested a maximum ca. 1695, we can 
estimate an average solar cycle length around the MM (from 
1636 to 1711) to be 9.5 + 0.5 years. If however, we assume that 
there were no additional solar cycle maxima around 1648 and 
1693, the average cycle length (max-to-max) would be 13.2+0.6 
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Fig. 3. Annual active day fraction for the three models. Error 
bars represent the 90% two-sided uncertainties. 




Years 


Fig. 4. Triennial active day fraction for the three models. Error 
bars represent the 90% two-sided uncertainties. Digital data is 
available in Table Q] 


years. However in this case, the length of individual cycles varies 
greatly, between 9 and 18 years. The estimated cycle length is 
similar to but so mew hat shorter than the results proposed by 
iMendozal d 19971) and llJsoskin et al.1 (1200 lh who suggested the 
cycle length of 10.5-11 years during the MM using sunspot ob¬ 
servations. Meanwhile, clustering of activity in » 20-year inter¬ 
vals (1650-1670,1670-1690, and 1690-1710) is also visible, in 
agreement with earlier results of the dominant 22-year period¬ 
icity during the MM (tlJsoskin et alJl200D . Note, however, that 
this clustering of activity could be also produced because of the 
scarcity of reliable data around 1648, 1669, and 1693. 

On the other hand, estimates of the cycle length based on 
cosmogenic 14 C data suggest much longer cycles during Grand 
minima (13-16 years). We note however that 14 C data cannot 
resolve individual cycles, because of the global carbon cycle 
attenuating high-frequency variability dRoth & Joosll2013l) . but 
rather yields the mean periodicity over the interval analyzed (e.g. 
Miva hara et all 120041) . Th is seeming contradiction between the 
results obtained here (cf. lMendozal1l997t lUsoskin et al.ll200~Tl) 
and from 14 C data can be potentially reconciled in a view of 
the possible inversion of the cycle phase in the cosmic ray mod¬ 
ulation during the periods of very weak activity like the MM 
(lOwens et al.ll2012l) . Thus, one or two cycles can be lost in the 
14 C data, due to forward and then reverse phase shifts in the be¬ 
ginning and end of the Maunder minimum, leading to a seem¬ 
ingly extended cycles in l4 C data. 


3.3. Sunspot numbers 


On average, the fraction of active days observed during MM 
was low, below 0.4 in the triennial data (Fig. El for ML and 
MO models, except for the year 1666 (MO model) which is 
however based on a single observation, and reaching up to 
> 0.7 in the most conservative MS model. We note that, for 
the normal cycles, the active day fraction is abo ut 1 00% ex¬ 
cept for the years around so lar minimum (Kovaltsov et al. 2004; 
lYaauero et alJl2012l 120141) . The value of Fa was never below 
0.15 for annual and 0.29 for triennial (see Fig. [5} during the 
period 1850-1995. Accordingly, such low values Fa even for 
the peaks during the MM correspond to (or are lower than) the 
minimum state of modern solar cycles. Therefore, although a 
cyclic activity during the MM is clear, at least during the core 
period, the sunspot cycles were weak, with the maxima being 
comparable to the modern cycle minima. We note that high so¬ 
lar cycle s of the magnitude 40-100 i n sunspot number as pro¬ 
posed bv lZolotova & Ponvavinl( ?015l) would unavoidably imply 
« 100% active day fraction (Vaquero e t al J|20 l~4l) during most of 
the years, which contradicts with the data (cf. Fig. [2}. 

In order to assess the sunspot number R from the acti ve day 
fraction Fa , we apply a method adopted from dKovaltsov et al.1 
l2004t Vaquero et al.l2012 , 20141) . For the annual data the relation 
was dKovaltsov et al. 2004): R = 19 • F^ 25 for F A < 0.5 and 
R = 2.1-exp (2.69 • F A ) for 0.5 < Fa < 0.8. The relation between 
triennial values R and F A is shown in Figure [5] for the period 
1850-1995. One can see that the relation is quite good for Fa < 
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exceeding 5 for the deep MM and 7 ca. 1705 (except for the 
very uncertain period ca. 1666 with the lack of observations). 
The most conservative MS model yields sunspot cycles below 
10 during the core MM and a possible relatively high cycle in the 
1650s, which is based on the lack of overlapping records from 
different observers, and about 15 ca. 1705. Anyway, even these 
very strict model suggests that the cycles were lower than 15-20 
in sunspot numbers, which is much lower than the present cycle 
#24 and an order of m agn itude lower than the very high cycles 
proposed by [Zolotova & Ponvavinl ([20151) . Considering the se¬ 
vere reduction of the statistics and a possible strong bias towards 
active days in the MS model, we believe it is not indicative for 
the true solar activity evolution during the MM and may repre¬ 
sent only the uppermost upper (maximum maximorum) bound. 


Fig. 5. Relation between triennial sunspot numbers and active 
day fraction for the period 1850-1995 using the Group Sunspot 
Number dHovt & Schattenlll998al) . The red curve is the best fit 
relation R = 33.6 • F l A 12 . 



Fig. 6. Annual (panel A) and triennial (panel B) sunspot num¬ 
bers reconstructed in the three models as denoted in the legends. 
Years with low statistics (Nr < 10) are not shown. Stars indi¬ 
cate that the sunspot number cannot be assessed from the ac¬ 
tive day fraction (see text) and is greater than 18/23 for the an¬ 
nual/triennial data. 


4. Conclusions 

Using three models of different level of conservatism to treat 
generic ”no-sunspot” statements, we have created a database 
of reliable sunspot observation around the Maunder minimum 
(1637-1715) and revised the sunspot cyclic activity over that pe¬ 
riod. We show that: 

1. A large number of no-spot records, corresponding to the so¬ 
lar meridian observations, may be unreliable in the HS98 
database. 

2. The active day fraction remained low (below 0.3-0.4) 
throughout the MM, indicating the low level of sunspot ac¬ 
tivity. 

3. The solar cycle appears clearly during the core MM with 
maxima at 1657, 1675, 1684, 1705 and possibly 1666. 

4. The length of the solar cycle during the MM appears shorter 
(9+1 years) in comparison with the standard 11-year solar 
cycle, but there is an uncertainty in that. A ~ 20-year clus¬ 
tering of activity is also observed. 

5. The magnitude of the sunspot cycle during MM is assessed 
to be below 5 (10 in the most conservative model) in sunspot 
numbers. The exact level is hardly possible to determine but 
it is below 10. 

6. Hi gh solar cycles during the Ma under minimum, as proposed 
by [Zolotova & Ponvavinl (12015i) . contradict with the data. 

We note that this is an uppermost upper (maximum maximo¬ 
rum) bound for solar activity during MM because of a possible 
selection bias (particularly important in the MS model), and the 
true level of activity may be smaller than that. 

In any case, only a thorough review of each record and each 
solar observation during the MM can make it possible to reveal 
the best picture of solar activity during this period. Therefore, we 
encourage researchers (especially Latin scholars) to query and 
analyze the old texts to understand how the observations were 
made and the true level of solar activity they indicate. 


0.8 (with the only outlier related to the period 1954-1956 which 
corresponded to the growth phase of the highest solar cycle #19) 
and can be well approximated by a dependence R - 33.6 • F X A 12 . 
The relations loosens for Fa > 0.8 and is lost completely with 
the active day fraction approaching unity. Thus, the active day 
fraction is a good index of sunspot activity until it reaches 0.8. 

Using these dependencies we have evaluated the sunspot 
numbers during the period analyzed, as shown in Figure [6] One 
can see that the sunspot numbers appear below 2 during the deep 
MM (1645-1700) and 7 ca. 1705 in the least conservative model 
ML. The optimum MO model yields the sunspot number not 
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Table 1. Triennial statistics of sunspot day occurrence for the three models considered here (see text for definition). Columns are: #1 
- central year of the triennial interval; #2 - number of active days N 4 within the interval; #3, 7 and 11 - number of total observational 
days Nr considered in the three models, respectively; #4, 8 and 12 - lower 90% bound of the active day fraction, for the tree models, 
respectively; #5, 9 and 13 - median active day fraction, for the tree models, respectively; #6, 10 and 14 - upper 90% bound of the 
active day fraction, for the tree models, respectively. 


Year 

n 4 

N r 

Flow 

ML 

Fmed 

F up 

N r 

Flow 

MO 

Fmed 

F U p 

N r 

Flow 

MS 

Fmed 

F up 

1639 

9 

10 

0.636 

0.851 

0.965 

10 

0.636 

0.851 

0.965 

9 

0.742 

0.932 

0.994 

1642 

135 

137 

0.956 

0.980 

0.992 

168 

0.752 

0.800 

0.843 

106 

0.956 

0.980 

0.992 

1645 

99 

171 

0.521 

0.576 

0.633 

167 

0.533 

0.590 

0.647 

99 

0.972 

0.992 

0.998 

1648 

0 

0 

N/A 

N/A 

N/A 

0 

N/A 

N/A 

N/A 

0 

N/A 

N/A 

N/A 

1651 

2 

20 

0.039 

0.123 

0.267 

20 

0.039 

0.123 

0.267 

3 

0.248 

0.613 

0.900 

1654 

31 

201 

0.120 

0.154 

0.195 

355 

0.068 

0.087 

0.109 

33 

0.827 

0.921 

0.974 

1657 

17 

170 

0.070 

0.100 

0.140 

52 

0.234 

0.329 

0.437 

27 

0.471 

0.622 

0.759 

1660 

66 

499 

0.114 

0.131 

0.151 

501 

0.114 

0.131 

0.150 

130 

0.439 

0.506 

0.573 

1663 

0 

168 

0 

0.003 

0.015 

0 

N/A 

N/A 

N/A 

4 

0.010 

0.128 

0.448 

1666 

1 

320 

0.001 

0.003 

0.011 

1 

0.223 

0.706 

0.974 

48 

0.006 

0.032 

0.089 

1669 

0 

541 

0 

0 

0.003 

0 

N/A 

N/A 

N/A 

109 

0 

0.005 

0.024 

1672 

47 

625 

0.064 

0.074 

0.086 

459 

0.086 

0.101 

0.121 

262 

0.148 

0.179 

0.215 

1675 

47 

561 

0.070 

0.082 

0.097 

196 

0.197 

0.239 

0.287 

108 

0.363 

0.434 

0.509 

1678 

30 

538 

0.045 

0.055 

0.068 

415 

0.057 

0.071 

0.089 

121 

0.193 

0.248 

0.312 

1681 

30 

516 

0.047 

0.057 

0.070 

738 

0.034 

0.039 

0.047 

146 

0.159 

0.205 

0.260 

1684 

72 

683 

0.093 

0.104 

0.116 

264 

0.235 

0.272 

0.312 

186 

0.334 

0.386 

0.440 

1687 

42 

697 

0.051 

0.058 

0.068 

432 

0.080 

0.096 

0.116 

190 

0.179 

0.221 

0.268 

1690 

12 

839 

0.011 

0.013 

0.016 

362 

0.022 

0.032 

0.047 

374 
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